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Abstract. We present observations of the formation of plate- like snow crystals that provide evi- 
dence for an edge-enhancing crystal growth instability. This instability arises when the condensation 
coefficient describing the growth of an ice prism facet increases as the width of the facet becomes 
narrower. Coupled with the effects of particle diffusion, this phenomenon causes thin plate-like 
crystals to develop from thicker prisms, sharpening the edges of the plates to micron or sub-micron 
dimensions as they grow. This sharpening effect is largely responsible for the formation of thin 
plate-like ice crystals from water vapor near -15 C, which is a dominant feature in the snow crystal 
morphology diagram. Other faceted crystalline materials may exhibit similar morphological growth 
instabilities that promote the diffusion-limited growth of plate-like or needle-like structures. 

1 Introduction 

The formation of complex structures during solidification often results from a subtle interplay of 
non-equilibrium, nonlinear processes, for which seemingly small changes in molecular dynamics at 
the nanoscale can produce large morphological changes at all scales. One popular example of this 
phenomenon is the formation of snow crystals, which are ice crystals that grow from water vapor in 
an inert background gas. Although this is a relatively simple physical system, snow crystals display 
a remarkable variety of columnar and plate-like forms, and much of the phenomenology of their 
growth remains poorly understood. (For a review of the physics of snow crystal growth, see [1].) 

Observations of snow crystal growth dating back to the 1930s [2l[T] reveal a complex and puzzling 
dependence on temperature. At pressures near one bar and water vapor supersaturation levels often 
found in clouds, for example, ice crystals typically grow into thin plate-like forms near -2 C, slender 
columns and needles near -5 C, very thin plates again near -15 C, and columns again below -30 C. 
The observed variations in snow crystal structure with temperature and supersaturation are often 
displayed in a snow crystal morphology diagram, as shown in Figure [1] Despite more than a half- 
century of study, we still do not understand the basic physical mechanisms that are responsible for 
the unusual temperature-dependent morphologies of growing ice crystals [U [3j |4] . 

Our focus here will be on the growth of thin plates near -15 C, as the largest and thinnest crystals 
form in this temperature region, making this a dominant feature in the morphology diagram. Aspect 
ratios of La/Lc ~ 100 or larger are commonly found with atmospheric ice crystals, where La is the 
crystal size along the a-axis of the crystal (roughly equal to the diameter of the plate), and Lc is 
the size of the plate in the c-direction (equal to the plate thickness). Large stellar plates typically 
have plate diameters of 2-4 mm with thicknesses of 10-20 /j,m [2]. 

In modeling the growth of ice crystals from water vapor, we write the perpendicular growth 
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Figure 1; A snow crystal morphology diagram showing the variations in growth forms as a function 
of temperature and supersaturation for crystals grown in air at a pressure of one bar (from [1]). 
The water saturation line shows the supersaturation of supercooled liquid water with respect to ice. 
Note that thin plate-like crystals appear only in narrow temperature ranges near -2 C and -15 C. 



velocity of a surface in the usual Hertz-Knudsen form 

^{Psurf - Psat) 
V = a r=^= 

v27rmfcT 



(1) 



where Psurf is the water vapor pressure at the growing surface, Psat is the water vapor pressure in 
equilibrium with ice, m is the mass of a water molecule, and f2 is the molecular volume [T]. This 
relation defines the condensation coefficient, a, that parameterizes the attachment kinetics at the 
ice surface, as well as the kinetic velocity Vkin, where cTsurf = (Psurf — Psat) I Psat is the water vapor 
supersaturation at the surface. In general a will be different for the two principal facets of ice, and 
for fast kinetics on non- faceted surfaces we expect a — >■ 1. 

If the supersaturation had a constant value independent of location near the ice surface, then 
the growth of a thin plate would imply aprism. / otbasai ~ Vprisva/vbasai ~ La/L^ > 1. Diffusion 
modeling shows that the supersaturation at the edges of a thin plate is not substantially higher than 
at the middle of the basal facets, and for a large fast-growing plate the minimum supersaturation is 
actually found at the edges of the plate [1] . Thus diffusion-limited growth will generally not promote 
the growth of thin plates in the absence of a large aprism/ctbasai ratio, and in many cases diffusion 
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requires that we must have aprism / oihasai > La/Lc, which further accentuates the difference in 
attachment kinetics on the two facet surfaces. 

The problem with this simple inference is that measurements of ice crystal growth rates for 
faceted surfaces have yielded values of aprism /ctbasai that are far smaller than what is needed to 
form thin plate-like crystals. For example, numerous researchers have found that ice crystals grow 
into roughly isometric shapes at low background pressures, where effects from particle diffusion are 
reduced [5l[6j[7], suggesting that aprism / otbasai is near unity for the conditions in these experiments. 
Since the growth of thin plates requires aprism / a^asai ^ 1, the discrepancy between low-pressure 
and high-pressure measurements necessarily requires a nontrivial explanation. 

In an earlier paper we proposed a mechanism to solve this problem that we called structure- 
dependent attachment kinetics (abbreviated SDAK) [5|. In essence, we proposed that aprism depends 
strongly on the structure of the prism facet surface, such that aprism increases when the width of 
the prism facet decreases, becoming near unity when the width approaches atomic dimensions. For 
the case of a thin plate, the width of the facet surface on the edge of the plate is approximately 
w w \/ aR, where a is the molecular step height and R is the effective radius of curvature of the 
edge of the growing plate [S]. Taking i? « 0.5 /xm for a typical thin plate yields w « 40a in the case 
of ice. We proposed that the SDAK mechanism leads to an edge-enhancing growth instability that 
promotes the formation of thin plates. 

Below we present observations of a hysteresis behavior in the growth of thin plates that strongly 
supports our proposed SDAK mechanism. Quantitative measurements reveal that aprism is indeed 
substantially larger on the edge of a thin plate than on a flat prism facet. This behavior then leads 
to an edge-enhancing growth instability that is responsible for the formation thin plate-like crystals 
near -15 C. 

2 Observations 

2.1 Intrinsic Growth Rates of Ice Facets 

We recently published improved measurements of the growth rates of faceted basal surfaces as 
a function of temperature [9 , indicating atasai ~ &^p{—o'basai/<^surf) over a broad temperature 
range, and specifically abasai = Cfcasa/,o ~ 2.0 percent at -15 C. Similar measurements of prism 
growth rates (lOj give the same functional form for aprism at -15 C, and for flat prism facets we 
have fjprism = o'prismfi ~ 4.2 percent. These measurements were made in a background of air 
at a pressure of approximately 0.03 bar, since lower pressures reduce effects from diffusion-limited 
growth. Similar measurements at 1 bar yield a values that are consistent with the low-pressure data 
[9]. These measurements are roughly consistent with previous data, although several improvements 
in our experimental techniques have yielded (Jprism,o values that are higher than quoted in pj. 

These data indicate that aprism / abasai ~ exp(— Acr/cr) with Act = 2.2 percent at -15 C. Since 
this ratio is less than unity for all supersaturations, it suggests (in the absence of the SDAK mech- 
anism) that the intrinsic growth form at -15 C should be a columnar crystal. Demonstrating this 
intrinsic columnar growth directly in low-pressure observations is difficult, however, since any ice 
surface contacting a substrate may be affected by substrate interactions. Generally we have found 
that these interactions reduce the nucleation barrier and yield untrustworthy growth measurements, 
especially at low supersaturations when the intrinsic 2D nucleation rate is low. Nevertheless, when 
one basal facet is flat against the substrate (as in j9j), we do occasionally observe the growth of tall 
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columns. Since substrate interactions cannot reduce the nucleation barrier all along a tall column, 
our interpretation of these observations is that sometimes the substrate interactions are greatly re- 
duced, due to some impurity coating on the substrate that we have not yet determined. When this 
happens we witness the intrinsic growth rates on all facets, giving columnar crystals. More typically 
substrate interactions increase the growth of the prism facets via enhanced nucleation where the 
crystal contacts the substrate, resulting in more isometric ice prisms. 

Our overall conclusion from this series of measurements is that the intrinsic attachment coeffi- 
cients for large facet surfaces are adequately described by the functional form above with (Tbasai,o ~ 
2.0 percent and Oprismfl ~ 4.2 percent at -15 C. Growth at low background pressures yields mainly 
large faceted surfaces, and the available evidence suggests that the intrinsic growth form at -15 
C may indeed be a columnar crystal. Additional low-pressure measurements, preferably with free 
falling or levitated crystals [11] to avoid substrate interactions, are needed to confirm this surpris- 
ing conclusion. The fact that thin plates are usually observed at -15 C is then entirely because of 
the SDAK mechanism. The resulting edge-enhancing instability appears only at higher background 
pressures, since the instability depends in part on the effects of diffusion-limited growth. We now 
examine additional measurements supporting this claim. 

2.2 Plate-on-Pedestal Growth 

We begin with a detailed description of the growth of a representative crystal at -15 C in air at a 
pressure of one bar. Using the apparatus described in [12], we began the experiment by dropping a 
small hexagonal plate crystal onto a temperature-controlled sapphire substrate. The initial radius 
of this crystal (taken to be half the distance from one prism facet to the opposite facet) was 13.5 
/im, and the initial thickness was 2.5 /im. The radius was determined by direct optical microscopy 
while the thickness was measured by broad-band interferometry, as described in [12 . Initially the 
temperature T/^ of an ice reservoir above the crystal was equal to the substrate temperature Tsubst , 
so the supersaturation was equal to zero. We then slowly increased the ice reservoir temperature 
and thus the supersaturation above the test crystal, causing it to grow. 

The crystal radius, its thickness, and the growth velocities as a function of time are shown in 
Figure [U along with the supersaturation above the crystal CToo • Note that cr oo was changed with 
time during the experiment by changing Tm, giving croo(t) as shown in the Figure, and this was the 
only external experimental parameter that was changed during the run. Since the distance from the 
substrate to the ice reservoir was much greater than the size of the test crystal, a^o can be assumed 
(for modeling purposes) to be the supersaturation at a hemispherical boundary far away from the 
crystal. 

As (Too was increased, the ice prism initially maintained its simple shape, namely that of a thick 
hexagonal plate. At roughly t w 160 seconds, the morphology of the crystal changed as a thin plate- 
like crystal began to grow out from the upper edge of the prism. We call this "plate-on-pedestal" 
growth because the thin plate grew outward above the substrate, supported by a central pedestal, 
like one half of a capped column crystal [1^ . 

Figure [3] shows images of the plate-on-pedestal crystal at various stages during its growth. The 
morphology of the thin plate during its rapid growth phase is essentially that of a sectored plate 
crystal jlj , and from additional observations we determined that the ridges along the a-axes of the 
crystal were on the underside of the plate, on the basal surface nearest the substrate. The top of the 
plate appeared to be a flat basal facet surface without ridges. In Figure [2] the radius refers to that 
of the top plate, while the thickness gives the overall thickness of the crystal from the substrate to 
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Figure 2: An example of "plate-on-pedestal" ice crystal growth. Top: The radius and thickness of 
the test crystal as a function of time. One basal surface of the crystal was resting on a substrate, 
and the test chamber was filled with air at a pressure of one bar. Middle: The basal Vh and prism Vr 
growth velocities. Bottom: The supersaturation far above the crystal CTqo and the supersaturation 
at the basal surface Csur/- The latter quantity was determined from the basal growth velocity using 



Equation[Tl assuming the measured condensation coefficient a\,asai 
in the text. 



exp(-crf,asai,o/o's«r/) described 
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Figure 3: An example of "plate-on-pedestal" ice crystal growth. Each image is labeled with the time 
in seconds, and the crystal sizes and growth velocities at each time are given in Figure [2j Between 
t — 130 and t — 180 a thin plate began growing from the top edge of the ice prism. Throughout its 
subsequent growth this plate did not contact the substrate except via its supporting pedestal. By 
i = 301 the outer edge of the plate had begun to thicken as the radial velocity diminished. 
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Figure 4: A plot of the basal growth velocity Vh and prism growth velocity Vr of the test crystal 
(see also Figure [2]). By t = 166 the thin plate had formed and the prism velocity quickly increased. 
After t = 178 the supersaturation was reduced and the edge of the plate thickened as both Vh and 
Vr became smaller. 



the top basal surface. Once CTqo was reduced, the edge of the thin plate thickened, and eventually a 
thick "rib" appeared around the circumference of the plate, as seen in the last image in Figure [31 

It is instructive to examine a Vh-Vr plot of the crystal growth shown in Figure HI as this exhibits 
a characteristic hysteresis behavior. In the early stages of the experiment, both Vh and Vr increased 
slowly as aco was increased, in this case with Vr ~ 1.5vh- Since the prism facets were in contact with 
the substrate during this phase of the growth, Vr was likely influenced by substrate interactions that 
increased nucleation on the prism facets and thus increased Vr compared to that for a free-standing 
prism facet surface. Had substrate interactions been absent, we would expect Vr < Vh because 
o'prismfl > cTbasaifl- By t = 166 the thin plate had begun to grow out from the top edge of the ice 
prism, which changed the structure of the prism surface from a broad facet surface to a sharp edge. 
According to our interpretation of the subsequent growth behavior, the sharp edge of the prism 
surface resulted in a higher aprism (and thus a higher Vr) because of the SDAK effect. This brought 
into play a growth instability that sharpened the edge of the plate further, leading to a rapid increase 
in Vr while Vh (and thus the supersaturation near the crystal surface) remained essentially constant. 
Around t = 178 the supersaturation aoo was reduced as shown in Figure [21 which caused Vh to drop. 
Because the edge of the plate was still sharp, however, aprism (and thus Vr) remained high. As aoo 
was reduced further, eventually the edge of the plate thickened and aprism dropped back down as 



7 




100 200 300 400 500 600 700 
time (seconds) 



Figure 5: Estimates of a surf and aprism as the test crystal grew, using the basal facet as a "witness 
surface" as described in the text. 

it approached the intrinsic value for a flat facet. 

The supersaturation field around the growing crystal was determined by diffusion, and we had no 
direct measure of isurf- However we can use the flat basal facet as a "witness surface" to estimate 
Osurf as follows. Since we know ai,asai{(^surf) — G'>^Pi—<^basai.o/<^surf) from other measurements of 
the growth of large basal facets, we can determine Usurf from this and the measured basal velocity 
Vh using Equation [T] The result is plotted in the bottom panel of Figure [2J We then further 
assume that Uprism — (Aprism I surf) and use Vr and Equation [T] to estimate aprism, with the 

result shown in Figure [5] Although somewhat crude, this quantitative estimate clearly confirms our 
hypothesis that the rapid increase in Vr when the thin plate formed is indicative of a rapid increase 
in aprism- Modeling the attachment coefficient as aprism — e:>q3{— aprism /a surf), this indicates that 
aprism <C aprism,o ou the edge of the thin plate, as expected from the SDAK mechanism. 

A key point in this analysis is an examination of what happened to this test crystal around 
t =^ 166. What caused Vr to suddenly increase while Vh remained essentially constant, as seen in 
Figure S]? Essentially the only change that occurred was that a plate- like appendage began growing 
on the top edge of the prism facet, as seen in Figure [31 This alone would not substantially alter 
the diffusion field around this small crystal, as is evidenced by the fact that remains roughly 
unchanged. Since the diffusion field did not change abruptly near t — 166, the rapid increase in Vr 
means that aprism increased threefold while asurf remained essentially constant. Our conclusion, 
therefore, is that aprism depends on the structure of the prism facet, which confirms our original 
hypothesis of structure-dependent attachment kinetics. 
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Figure 6: A second example crystal grown at -15 C, showing a more gradual increase in CToo with 
time during the early stages of growth. 
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Figure 7: Top: A plot ol the basal growth velocity Vh and prism growth velocity Vr of the test crystal 
shown in Figure [6l Bottom: Three images of the crystal at different times during its growth. 
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Figures [S] and [7] show another example where (Too was increased more gradually as the plate 
developed. In this case the plate initially formed at a lower Vh compared to the previous example, 
and the maximum Vr was lower. Again the Vh-Vr plot shows significant hysteresis, indicative of the 
SDAK mechanism. Here we see that aprism increased by a factor of four as the thin plate developed, 
while during this same time Vh (and thus a surf) declined slightly. Modeling a prism as before yields 
a broad dip with a minimum value of aprism ~ 0.8 at t w 450 seconds. 

We observed similar plate-on-pedestal growth behaviors at other temperatures ranging from -12 
C to -17 C. The hysteresis in the Vh-Vr plots was most pronounced at -15 C, and we observed that 
the minimum Vh needed to produce a transition to plate-like growth was approximately two times 
higher at -17 C than at the other temperatures. Overall these observations are consistent with what 
we know from the morphology diagram, namely that the growth of thin, plate-like crystals is most 
pronounced near -15 C, as shown in Figure [TJ 

3 Modeling 

It is also instructive to examine computer models of growing crystals as aprism and other parameters 
are varied, in order to compare with our observations. To this end we used a cylindrically symmetric 
cellular automata model [13] for which the functions abasaiiasurf) and aprism{a sur f) were inputs in 
the code. For these computations we used abasai[asurf ) = exp{-abasai,o/asurf ) with abasai,o = 2.0 
percent and aprism{asurf) — exp{— aprism /asurf) with a constant aprism for each model. These 
models cannot reproduce the full SDAK growth instability, or the observed hysteresis, because 
aprism was assumcd to have a constant value independent of the structure of the crystal. Therefore 
we used the models mainly to examine the transition to plate-like growth as a function of aprism 
and (Too- 

Figure [5] shows one example of a set of models as aprism was varied. The figure shows that the 
models produced thin plate-like crystals growing from the initial pedestal only when aprism was quite 
low, confirming the low aprism values we inferred from our measurements. From this and additional 
model calculations we reached two conclusions: 1) both the observed plate-on-pedestal morphology 
and the measured Vr values were approximately reproduced with aprism values that were consistent 
with those inferred from our experimental data, and 2) the initial growth phase with no thin plate 
was reproduced only with larger aprism values, again roughly consistent with our measured values. 
Although these constant-CTprism models cannot reproduce the full behavior of our observed crystals, 
they nevertheless support our overall conclusions regarding the SDAK mechanism and the resulting 
edge-enhancing growth instability. 

To produce a full model of the SDAK instability, it would be necessary for aprism to vary with 
the structure of the crystal, as this is an inherent feature of structure-dependent attachment kinetics. 
One possibility would be to let aprism{a sur f) — i^pi—aprismiR) / asurf) where R is the radius of 
curvature of the edge of the crystal plate. The model would assume aprism{R = oo) = aprism.Oi where 
aprismfi is the measured value for flat prism facets. At smaller R, aprism would become smaller, 
as seen in Figure [51 Realizing such a dynamical model would require some guesses for the correct 
functional form for aprism{R), as this function is not determined from our measurements or from 
theoretical considerations. Without more constraints on this and other model inputs, we believe a 
full dynamical model would be of limited use at this time, other than as a proof-of-principle that the 
SDAK instability can exhibit a hysteresis behavior that approximates what we observed. However, 
additional experiments like those described here may yield a useful measurement of aprism{R), so 
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Figure 8: Models of plate-on-pedestal growth as aprism is varied. Shown are cross-sections of 
the cylindrically symmetric crystals after 20 seconds of growth. Brightness around the crystal is 
proportional to supersaturation. For each of these models the initial crystal radius and diameter 
were both 5 iim, CToo = 5 percent, abasai = e^pi-crbasai,o/crsurf), and aprism = ewi- (Aprism /(^ surf)- 
Each model has the same scale, and the (Jprism = plate has a final radius of 26 /xm. The prism 
growth velocities for the final crystals shown were Vr = (250, 300, 450, 600, 650, 900) nm/sec, from 
top to bottom. 
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additional progress in quantitative dynamical modeling is certainly possible. 

4 Discussion 

To summarize, the data and modeling presented above provide strong evidence for structure- 
dependent attachment kinetics and an edge-enhancing growth instability at -15 C. Our ice crystal 
growth measurements and modeling indicate that: 

1) For the growth of flat facet surfaces near -15 C, the condensation coefficients for both the 
basal and prism facets are well described by the functional form a{a) = exp(— cr/acet.o/o"), where a 
is the supersaturation at the surface and fJfacetfi is a constant. 

2) For flat facet surfaces we have measured atasaifi ~ 2.0 percent and aprism,o ~ 4.2 percent at 
-15 C. Because abasai.o < ctprism.o, this suggests that the intrinsic growth morphology at -15 C is a 
columnar prism. 

3) The measurements suggest that the values of abasai,o and aprism,o are independent of back- 
ground air pressure, at least for the range 0.01 - 1 bar. We also expect there would be no intrinsic 
pressure dependence in the attachment kinetics from theoretical considerations, since these pressures 
are too low to substantially affect the ice surface. 

4) While chemical contamination of the ice surface is always present at some level in experiments, 
we have found that contamination levels must be quite high to affect ice growth [14j . Furthermore, 
chemical contamination typically suppresses the growth of thin plates at -15 C [14], suggesting that 
contamination is not an important factor in the observations described here. 

5) Our observations of plate-on-pedestal growth show a rather abrupt transition to plate-like 
growth. Using the basal facet as a "witness surface" allowed us to estimate aprism as described 
above. During plate growth we found aprism ^ ^prismfi- 

6) Computer modeling confirms that aprism ^ aprism,o is necessary for plate-like growth. 

7) Hysteresis in the Vh-Vr plots requires that aprism{a sur f) is not a simple single- valued function 
independent of crystal morphology. Instead, aprism must be large at the edge of a thin plate and 
smaller on the edges of thick plates or on flat facet surfaces. 

Together these considerations strongly support the SDAK mechanism we proposed earlier, and 
indicate the presence of an edge-enhancing growth instability. Beginning with a crystal prism con- 
taining only large facets, we initially have aprism ~ Cprism.o and abasai ~ a^asaifi- At low pressures 
the prism and basal facets remain large and the SDAK mechanism is absent, giving the roughly 
isometric growth that has been observed. But at higher pressures, diffusion causes the top edge of 
the crystal to grow a bit faster than other parts of the crystal, so the edge sharpens slightly. This 
sharper edge leads to a reduced aprism-, which causes the edge growth to increase further. A positive 
feedback results: as aprism reduces, the edge grows faster and becomes sharper, thus reducing aprism 
more, and so on. This positive feedback is the hallmark of a growth instability, and the end result 
is the formation of a thin plate. 

An edge-enhancing instability driven by the SDAK mechanism is a relatively simple physical 
explanation that nicely explains all the observations. The existence of a growth instability of this 
nature may also explain why ice crystal morphologies change so abruptly with temperature in the 
morphology diagram. If the instability tips only slightly in the direction of fast prism growth, then 
thin plates will form. But if the same instability tips toward fast basal growth, then hollow columns 
(which, like plates, also have thin edges) become the preferred morphology. Even small changes in 
the underlying nature of the instability can result in large changes in the final crystal morphologies. 
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We do not yet understand the molecular mechanisms that produce the SDAK effect, although 
structure-dependent surface premelting may be playing a role. We also do not yet understand the 
temperature dependence of this instability, and why it is especially prevalent at -15 C. Nevertheless, 
it appears that this instability is responsible for the formation of thin plates at -15 C, which is 
a dominant feature in the morphology diagram. Additional observations may yield more insights 
into the underlying molecular processes responsible for shaping the growth of snow crystals under 
different conditions. 
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